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The effect of the electron temperature on both the light absorption and the scattering by metal 
nanoparticles (MNs) with excitation of the surface plasmon electron vibrations are studied in the 
framework of the kinetic theory. The formulae for the electroconductivity and polarisability tensors 
are derived for finite temperatures of an electron gas. The electrical conductivity and the halfwidth of 
the surface plasmon resonance are studied in detail for a spherical MN. Depending on the size of MN, 
the efficiencies of the light absorption and scattering with temperature changing are investigated. 
It is found, in particular, that the absorption efficiency can not only increased with a temperature 
decreasing, but decreases as well. Derived formulas allow to calculate analytically the various optical 
and transport phenomena for MNs with any spheroidal shape embedded in any dielectric media. 
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I. INTRODUCTION 

When the metal nanoparticle is illuminated with a 
laser beam by frequency which coincides with the fre- 
quency of collective electron oscillations in the MN, 
the surface plasmon resonance (SPR) is excited. The 
frequency and the width of SPR depends on the 
size, morphology, spatial orientation of MN and on 
the dielectric environment.^'^ Resonance light scatterers 
are employed in wide range applications ranging from 
surface enhanced Raman scattering ^2, near-field scan- 
ning optical microscopies j^ to bio-chemistry imaging;^ 
surface enhanced fluorescence jS subwavelength optical 
waveguidesj2. cancer therapy^ etc. 

The temperature effect on the optical and transport 
properties of the metal nanoparticles is very important 
for pure and applied science of the nanoparticles.— As a 
MN absorbs laser energy in a thin-surface layer, rapid lo- 
cal heating can occur at the surface. Thus, for a detailed 
analysis of the laser-light absorption or reflection at MN 
surfaces, the properties of MNs at electron temperatures 
must be studied. The temperature dependence of the SP 
resonance is crucial due to the number recent applications 
of noble MNs in computer chips, ^° thermally assisted 
magnetic recording, ^^ thermal cancer treatment,— and 
catalysis and nanostructure growth^i^ The use of nano- 
objects as temperature sensors^ and thermometersi^ is 
promising due to their small sizes and short thermal re- 
laxation time. 

Previous calculations the effect of the temperature on 
the different properties of the MNs were performed in 
the zero-temperature limitii^iii or in the interval 15-1000 
K onlyi^ii^"— While the optical properties of both the 
bulk metals^^Ti^ and the MNs in the low-temperature 
regime is now well studied, an understanding of the ef- 
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fccts of electron temperature (typically of order lO^i^T) 
on the plasmon modes has remained a challenge and is 
not as well understood as at the T — Q. Such an out- 
of-equilibrium situation can easily be induced in metallic 
nanoparticles^^ using ultrashort optical pulses. At time 
t — Q, only the electrons are excited by the laser pulse; 
then the excited electron would decay through light radi- 
ation, phonon excitation, and/or electron-electron scat- 
tering. In the most commonly used approach, the two- 
temperature modelj^i^i^ which is assumed that both 
electronic and ionic degrees of freedom are in thermal 
equilibrium conditions individually, but each degree of 
freedom has its own temperature. The energy exchange 
between the electrons and the lattice in this model lead 
to a time-dependent electronic temperature Te(t). 

The effect of electron temperature was not studied in 
detail to date because a broad temperature interval re- 
quires a usage of materials with high thermal stability. In 
Ref. I22II27I the improved Drude model when the electron 
concentration and/or a plasma frequency are directly de- 
pendent on the electron temperature were used. Only few 
works have so far addressed this problem. ■^^''^^ 

In this work, we focus on optical properties of metal 
nanoparticles at relatively high electron temperatures, 
whereas the ion temperature, Ti , is supposed to be much 
smaller,^'^ Different behavior with temperature of opti- 
cal lines for MNs with different radii and shapes are a 
particular motivation for our present study. The screen- 
ing coming from the surrounding matrix is taken into 
account through a constant dielectric function, £„. 

To calculate the effect of the electron temperature on 
the optical and transport characteristics of MNs, we will 
use the kinetic theory method which accounts for the 
electron scattering on the MN boundary. It is shown that 
the efficiency of absorption with temperature strongly de- 
pends on the size of MN. It is predicted that the effect 
of temperature on the MNs with different radii accounts 
for most of the possible optical properties in this inter- 
esting temperature regime. We present our theoretical 



results on the temperature dependence of the conductiv- 
ity and polarizability tensors related to the MNs with the 
spheroidal shape. 

The rest of the paper is organized as follows. The 
Boltzmann equation approach to the problem is pre- 
sented in Sec. II. Section III contains the study of the 
conductivity tensor in the spheroidal MNs at finite tem- 
peratures. In Sec. IV, we consider the effect of tempera- 
ture on the polarizability of MN. Section V is devoted to 
study of the light absorption and scattering crossection 
by MNs. The discussion of the obtained results is avail- 
able in Sec. VI, and Sec. VII contains the conclusions. 



II. BOLTZMANN EQUATION 

To account for the temperature effect on the conductiv- 
ity and polarization properties of the MNs, we will apply 
the Boltzmann equations approach. The advantages of 
this approach is that the obtained results can be applied 
not only to MNs with a spherical shape, but to strongly 
anisotropic spheroidal (needle-like or disk-shaped) MNs. 
Thus, it permits one to study the effect of the parti- 
cle shape on the measured physical values. Second, the 
Boltzmann equation method enable us to investigate the 
MNs with sizes smaller than the electron mean free path 
I. But in the case of sizes less than I, it gives the same 
results as known from other approaches. 

There exists a lower limit of the applicability of Boltz- 
mann method in the small radius limit when the particle 
size is comparable to the de Broglie wavelength of the 
electron, and the quantization of the electron spectrum 
starts to play an essential rolcii^ Practically, it is around 
a radius of greater or less than 2 nm. 

Let us consider the single ellipsoidal metal MN with 
semiaxes a, b, c that is irradiated by an electromagnetic 
(EM) wave, whose electric field is given as 



connected to Eq by the relation^S 



E = Eg exp[i(kr — ivt)]. 



(1) 



Here Eq is the amplitude of an electric field of a pump 
laser, uj is its frequency, k is the wave vector, and r and 
t describe the spatial coordinates and time. 

We limit ourselves with the case of Rayleigh scatter- 
ing, when the electromagnetic wavelength (from a pump 
laser) A ~ c/w is much larger than the diameter of the 
nanoparticle d {— max{a,b,c}). Then the electromag- 
netic field around the MN can be considered as homoge- 
neous. Placing the coordinate origin in the center of the 
particle, the above-mentioned assumption is written as 



fcr < 1. 



(2) 



The inequality ([5]) implies that the E field of the elec- 
tromagnetic wave can be considered to be spatially uni- 
form on scales of the order of a particle size such that all 
the conduction electrons move in-phase producing only 
dipole-type oscillations. The amplitude of such a field is 



E^'\0,^)/E, 



U) 



(w) = l + L,[e(w)/e„,-I], (3) 



where e{uj) is the dielectric permittivity of the MN, e^n 
is the dielectric constant of the adjacent medium (i.e., 
solvent), and Lj are depolarization factors in the j-th 
direction (in the principal axes of an ellipsoid) . 

The electric field Ein gives rise to high-frequency cur- 
rent inside the MN. To obtain the average density of this 
current over the MN, it is necessary first of all to cal- 
culate the electron velocity distribution function. The 
field Ein has an effect on the equilibrium electron veloc- 
ity distribution and thus determines the appearance of a 
nonequilibrium addition /i(r, v,t) to the Fermi distribu- 
tion function 



foie) = 



1 



cxp 






1 



Here fJ-{T) is the chemical potential at finite temperature 
T, fcs is Boltzmann's constant, e — ■mv'^/2 denotes the 
kinetic energy of an electron, u = | v| refers to the electron 
velocity, and m is the electron mass. As is well known,— 
the equilibrium function /o(e) does not give any input to 
the current. Accounting for both the time dependence 
of Eq. ([1]) and the inequality ©, the total distribution 
function of electrons, can be represented as 

/(r,v,t)-/o(e) + /i(r,v,t) = /o(e) + /i(r,v)e'"*. (4) 

The function /i(r, v) we will seek as a solution of the 
linearized Boltzmann's equation 

(.-zc.)A(r,v)+v^%::l + eE.„v^.O, (5) 



dr 



de 



where e is the electron charge. In Eq. ([5]) we have as- 
sumed that the collision integral 



(dfi/dt). 



-/i/r 



is evaluated in the relaxation time approximation (t = 
l/i/, 1/ refers to electron collision frequency). Strictly 
speaking, i^ = J^(r) is temperature dependent value 
which can be presented as 



HT) 



m Q 



e/T 



(e^ - 1)- 



■ dz, (6) 



where n,, is the electron concentration, Q is the Debye 
temperature, and K combines together factors depend- 
ing on the details of the Fermi surface geometry and scat- 
tering matrix elements. 

For simplicity, we assume as well that the vortex elec- 
tric field, induced by magnetic component of the external 
EM field, gives comparatively small input at the plasmon 
resonance frequencies and can be neglected in Eq. ([S]). 

What is more, the function /i(r, v) ought to satisfy 
the boundary conditions as well. These conditions may 



be chosen from the character of electrons reflection from 
the inner walls of the MN. We will adopt, as is usually 
done, the assumption of diffusive electron scattering by 
the boundary of MN. Then the boundary conditions can 
be presented in the form 



/(r, 



vjis 



0, 



<o. 



(7) 



where v„ is the velocity of the component normal to the 
particle surface. 

Along with diffusive scattering, the mirror boundary 
conditions at the nanoparticle surface were examined in 
the literature for electron scattering (see, e.g., Ref. \iOi ). 
In this case, each electron is reflected from the surface 
at the same angle at which it falls to the surface. In dif- 
fuse reflection, the electron is reflected from the surface 
at any angle. In order for the mirror mechanism to be 
dominant, the surface must be perfectly smooth in the 
atomic scale, since the degree of reflectivity of the bound- 
ary essentially depends on its smoothness. Practically, 
for a nonplanar border such smoothness is extremely dif- 
ficult to achieve. As was shown^^ the mirror boundary 
conditions give a small correction to the results obtained 
with the account of only the diffusive electron reflections. 
Therefore, we chose more realistic boundary conditions 
given by Eq. ([7]). 

The boundary conditions ([7]) in the case of an ellip- 
soidal MN, generally speaking, depend on angles, what 
complicates the solution of Eq. ^. 

It is rather easy to solve Eq. ([5|) and to satisfy the 
boundary conditions of Eq. (O if one passes to the trans- 
formed coordinate system, where an ellipsoid with semi- 
axes a, b, c (along the x, y, and z directions, respectively) 
transforms into a sphere of radius R with the same vol- 
ume: 



Xj = -^x'j, R={ahcY'^, 



(8) 



with i = 1,2,3, and di = a, d2 = b, d^ = c. A similar 
transformation should be made for the electron velocities: 






(9) 



Equation ([5]) and the boundary conditions ([7]) in trans- 
formed coordinate and velocity systems, can be rewritten 
as 



i.-^u)Mr'y)+V^■f'^''^^'^ ' - -'^^'^'^ 



eE.;„,v 



r oe 

/(r',v')|.'=fl = 0, rV<0. 



= 0, 
(10) 

(11) 



The first condition calls for zero equality of the electron 
distribution function at the nanoparticle surface, and the 
second one, rV < 0, means that an electron motion 
occurs only inside MN and the electron leakage through 
the NP surface is absent. 



Eq. (ITOl) presents the partial differential equation with 
the boundary conditions ([TT|). Using method of charac- 
teristics, one can obtain the solution to this equation in 
the form 

/i(r ,v ,t) = -e— — vEi„ ■. , 

oe u — iLo 

(12) 
where the characteristic ic(r',v') can be presented as 



ic(r',v') 



(R2-r'2)v'2 + (r'v')^ 



(13) 

The radius vector R determines the starting position of 
an electron at the moment <c = 0. The characteristic 
curve of Eq. (IT^ depends only on the absolute value of 
R and does not depend on the direction of R. 

It is reasonable to point out here that, in spite of the 
electric field remains spatially uniform inside the MN, 
the distribution function Eq. (|12p depends all the same 
on the coordinates in consequence of the requirement to 
obey boundary conditions Eq. (jlip . Owing to this depen- 
dence, other physical parameters averaged with /i(r', v') 
start to depend on coordinates too. Because, the phys- 
ical sensing have only parameters averaged over whole 
MN volume V , it is necessary to fulfill the integration 
over all coordinates inside of MN. One can see this on 
the example just below. 

Using the solution ([T^. one can calculate the density 
of a high-frequency current induced by the EM wave in- 
side the MN. Performing the Fourier transformation of 
Eq. ([H]), we obtain 



2e / m 

J(") = y (2^ 



d^j 



dVv'/i(r',v',c.). 
(14) 



III. ELECTRIC-CONDUCTIVITY TENSOR 

Let us introduce the tensor of electric conductivity 
(^ap{<^) using the relationship 



13=1 



(15) 



Then, in accordance with both Eqs. p^ and (HH), the 
components of this tensor can be presented as 



acp[uj) 



^fZL^l^ /// ^3„/ /// ^3„/, 

V \2TTnJ 



d'r' / / / d'v' V, 



evp 



dfo i-e-(''-*'^)*<=('^''"'')^ 



de 



(16) 



Electric-conductivity tensor CTq^ = ct^q, is the symmetri 
cal second-rank tensor. Because the integrand in Eq. ([T~~ 



is an odd function for a nondiagonal tensor components 
a 7^ /?, and the integration is conducting over the whole 
velocity space [—00, cx)], that retain only diagonal com- 
ponents in this equation. 

If we introduce a rectangular coordinate system and 
chose the fundamental directions x, y, z to be coincident 
with the three principal axes of the ellipsoids then we get 



'xx 











^m 











fj. 



We denote with ip and 9 the azimuthal and polar angles 
with respect to the ellipsoid rotation axis z, respectively. 
Note here that only diagonal terms with v^ — vp = v are 
retained after integration over all angles. As one can see 
from Eq. p^ . the conductivity becomes scalar quantity 
in this approximation. 

In the general case of an ellipsoidal-shaped MN, the 
electric conductivity is the tensor quantity (Eq. (J16p ). 
Integrating it over all nanoparticle coordinates in a solid 
angle (Kl = sin 6 dip dO gives (see Appendix) 



The tensor of electric conductivity, as can see from 
Eq. (fTB)) . is the complex value: 



1 

V 



d\' 



-{u—iLo)t (r ,v') 



-^{u,v')- (22) 



(Tq/3(w) 



'aP 



(CO) 



'at3 



[U). 



(17) 



The surface effect on the conducting phenomenon is 
described in Eq. ([T6| by means of the characteristic 
ic(r', v'). It accounts for the restrictions imposed on the 
electron movement by a nanoparticle surfaces. As one 
can see from Ec^. (IT^ . the value of tc is of the order of 
t' ~ R/vp, where vp is the Fermi velocity. This implies 
that the value reciprocal to tc will correspond to the vi- 
bration frequency between the particle walls. Hence, the 
inequality vtc 3> 1 indicates that the electron collision 
frequency inside the MN bulk would significantly exceed 
the one for an electron collision with the surface of MN. 
If this inequality is satisfied, one can to direct tc ^ 00. 
Then the exponent in Eq. (J16p vanish and finally we ar- 
rive at a standard expression for an electric-conductivity, 
such as following ^!^ 



aiuj) = 



1 ^;i 

At: v — iuj^ 



with 



2 Ae'^m^Vp 



iirh? ■ 



(18) 



(19) 



To ensure that, it is necessary to take into account that 
the energy derivative of /o in the zero approximation in 
the small ratio of ksT /ep {ep is the Fermi energy) can 
be replaced by 



— ^-S{s-sp), 



(20) 



and pass in Eq. (J16p to the integration over v in the 
spherical coordinate system 

d^ u -)■ dip sinOdO v'^ dv, 
00 



with the use of the formula 



vU{v^~vl)dv=^. 



(21) 



The complex ^ function entering in Eq. (E^ . has the 
form 

^{L.,v')^HL,,v')^^(l + -)e-^, (23) 

q \ qj 

with 

4 2 4 2R 

$(cj, v') ^ + — , q = q[uj, v') = —{y-iuj), (24) 

A q q'^ v' 

and v' (= <^v) is a "deformed" electron velocityii with 
the "deformation" coefficient <ij = R/dj. The last sum- 
mand in Eq. ([23]) represents the oscillation part of the ^ 
function and the first one refers to its smooth part. 

Accounting of Eq. ([^^ and only diagonal components 
in Eq. ()16p . it remains the integrals over all electron ve- 
locities. To calculate them, we pass to the spherical co- 
ordinate system with a z axis directed along the rotation 
axis of the spheroid (as we have done it above). Then 
Eq. (fT6|) can be rewritten as 



crjj(w) = 



^th'J 



3em 



1 



27r 



2{2Tihf v~iuj 



dip I sin 9 d9 



v^dvv]L^\^{Lo,v'), (25) 



where Vj is the j'-th component of the electron velocity, 
3 =x,y,z, with 



v%, = v^ sin^ e ■ (™%" ^) 



respectively. 

It should be noted that the "deformed" electron ve- 
locity entering in ^ function can be expressed through 
the electron velocity in Cartesian coordinate system as 
follows 



/ / cos^ 69 sin io 



62 



<s^e 



sm 



(26) 



Let us suppose that the particle in the matrix is mod- 
eled as rotationally symmetric ellipsoid {a = b = R±, 
c = i?|| ) with the symmetry axis along the z direction. 
The components of an electron velocity parallel {vh) and 
perpendicular {v±) to the spheroid revolution axis 



uii = Vz = V cost 



VI = xlvl + vl 



vsmO (27) 



play an important role in this case, and the v' cease to 
depend on the angle ip. 

Let us puss from an integration in Eq. (j25l) over elec- 
tron velocities to the integration over electron energies 



1 /2e 



V dv = — — de 



3/2 



(28) 



and account for that the integrals 



cos^ f — sin ip — TT. 



(29) 



Then, using Eqs. ^5^, ([281), and ^9^, one get for main 
components of the complex electric conductivity tensor 



where Eq. (p6)) was used. Because the spheroid semiaxes 
R± and Rh can be easily expressed through the radius of 
sphere R of an equivalent volume 



-R I ^= R 



R± 

Rn 



1/3 



R\\=R 



R± 

Rn 



-2/3 



(34) 



the energy e' does not depend on the particle radius R 
but depends only on the spheroid axes ratio. 

If we limited ourselves here with the case of low tem- 
peratures, then Eq. ([20)1 can be applied, and we found 



/ *(c^, 4) S{e - ep) e'/' de = fil^^'i'{oj, fi'e,), 
Jo 



(35) 



where /io is the chemical potential at zero temperature. 
Substituting Eq. (1551) into Eq. (I5T|) and accounting that 
the electron concentration at T = can be presented as 



no 



(2m/io 



,3/2 



3TT^h^ ' 



(36) 



it is easy to check that Eq. (PT|) transforms to the form 
known from, e.g., Ref. l42l 



ayy = a± , 



with 



o'(|)(w) 



3e2 V2to 
2^~fi3- 



7r/2 



/cos 6*' sin 61 






s^9' 



f.(.,4)(-t).^/^.e 



(30) 



dd' 



(31) 



We denote with d'{= 9 — 7r/2) the angle between the di- 
rection of an electron velocity and an axis perpendicular 
to the spheroid rotation axes. Here and below, the upper 
(lower) symbol in the parentheses on the left-hand side of 
Eq. pip corresponds to the upper (lower) expression in 
the parentheses on the right-hand side of this equation. 
The ^ function in Eq. (|3T|) depends now on both the 
energy e' and the angle 9' because the parameter q [see 
Eqs. (p4)) and (|26)) ] for a spheroidal particle becomes de- 
pendent on the angle 9 and can be determined as 



qi9,e), 



(32) 



where R\\ and R± are the semiaxes of the spheroid. The 
electron energy e' entering into the '3/ function, becomes 
dependent on angle 9 too and for a spheroid can be pre- 
sented as 



o2fl' 



eR^ 



sin^ 9' 



Rl 



R 



(33) 



A. Conductivity of a spherical MN 

For particles with a spherical shape, the electric con- 
ductivity becomes a scalar quantity, and one can put 



Rl 



Ra 



R, e' = e in Eqs. 



and 



then q 



and the vjf function cease to dependent on the angle 9' 

7r/2 



cos^'sin^^'X 1 

icos3 0' 7 3' ^ ^ 



and Eq. (|3T|) reduces to the form 
e^V2m 



r^sphi^) 



2n'^h?{v - ii^) 



r*<-<-t)^""*' 



(38) 

Let us make the following change of variables in 
Eq. dSHl) 






dij 



Then 



dfo _ e-" 

de (l + e-'')2 



de 



drj, 



(39) 



(40) 



and the integral in Eq. (|55)) can be presented as 



*(c.,77)(M + 7yfcBT)3/2^_-_^dr? = /, (41) 



-M/(fcBT) 



where 



Using Eqs. (|39l)-(|48l), Eq. (|38l) transforms into 



^ ' " 3 R{v-iuj) 2R^{v-iujf 

R^iv-iLof \ 2R{v^iLo) ' 



within 



^ = Ai(r) - A^o 



1 



TT^ /fcsT 



12 V /io 



(42) 



(43) 



Since the input of the integrand tends to zero as 77 — )■ 
— c», the lower limit in this integral can be extended to 
the — cx). The derivative {—dfo/ds) has a maximum at 
the point e = /i. So, it is conveniently to expand the 
product of ^ • (/i + riksT)^'^ = S(r/) in the powers of 77 



E{rj) c, ^3/2^(0) 



JfcBTVA^*(0) + Ai'/2^'(0) 



-^ [3(A:bT)2*(0) + 12fcBrAi*'(0) + 4Ai2l'"(0)] r/^ 
OM', (44) 



with 



^{0) = *(c^,77)|^^o, *'(0) = ^'i'{co,v)U^o, 



*"(o) = a;^*(^'^)l';-o- 



(45) 



Note, that ^(0) at v' = vp coincides with the ^(w) de- 
fined by Eq. (PI) . 

Then, the integral (jH]) can be presented as the sum of 
integrals 



^3/2^(0) /i 



HksTT 
4V7^ 



-fcsrV7I*(0) + Ai3/2^'(0) 



*(0) + SkeT^^'iO) + ///2*"(0) 



/a, 
(46) 



where the symbols /i , I2 and /a represent the integrals 

00 00 

e"'' dr] _, j_ /" *? 6~'' d?7 



h = 



(l + e-'')2 



1, /2 = 



(1 + e--iy 



= 0, 



and 



/a = 



77 e '' (i?7 TT^ 



(l + e-'')2 



(47) 



(48) 



<7sphi^) 



Here 



e -y/2r7T,/i 
27r2?i3(j,„i^) 



vr^ /fcsT 



8 V M 



*(0) 



[3kBT^'{0) + fi^"{0)]+0[T]]H. (49) 



*'(0) = 






-l + 4--e--('l + l + % 



9m 



9m 3 



(50) 



*»^^ 



9^ 



4e" 



9m 



^^ 2 +2g,, + 2q2 



with 



9m = 



2R{iy - iuj) 



(51) 



(52) 



Substituting Eqs. dSHll-lISIl) into Eq. (gH]), we obtain 



asph{^,T) 



me^^R 



TT^h? 



It^ 



with 



ax 



[*(0)+aT*T(w)], (53) 



(54) 



^t{^) = $t(w) 



^2 


(ksT' 


]' 




6 


\ M / 


1 ' 




8 
9m 


24 

+ — 


1 + 
V 


1 
9m 



-9k 



$T(a;) = 1 



4 24 
9m 9^' 



, (55) 



(56) 



When the temperature T — > in Eq. ([M)) . then 
Eq. ([53)) is reduced to the form 



0'sp/i(w) 



2..^/2^/2^ *(0) 



e Mo 



27r2n3 



(57) 



known from other calculations (see, e.g., Ref. |44 

The conductivity of MN given by Eq. ((55)) is a complex 
value. Practically, we need to know both the real and the 
imaginary parts of it. 

Further analytical calculations are possible for some 
particular cases. In the case of a spherical MN there 
are three actual frequencies that are considered usually: 
the frequency of an incident electromagnetic field w, the 
collision frequency of electrons in the particle volume v, 
and the vibration frequency between the particle walls 



i/g — vp/i^R) (if the particle size is less than the elec- 
tron mean free path). When i' > v^, the mechanism of an 
electron scattering in the bulk is dominated, and an elec- 
tron scattering from the particle surface gives only small 
corrections of the order of i^s/i^. But we are interested 
in the case when the mechanism of the surface electron 
scattering dominates, which corresponds to j^ < i^g. 

Let us suppose that the mean free path of a conduction 
electron I is much greater than the particle size. In this 
case, an electron scattering occurs mainly from the inner 
surface of the MN. The electrons oscillate between the 
walls of the particle with the frequency I's- We restrict 
ourselves here only to the case, when v <^ Vs- In this 
limit case, to a first approximation, formally one can put 
j^ — > 0. If one introduce the designation 



qi 






(58) 



then for the real and the imaginary parts of the ratios 
"^ /{v — iuj) and $T/(i^ — iw), we obtain 



Re 

Im 
and 



* 



* 



iy->0 



1/^0 



1 

UJ 

4 

UJ 



2 4 4 

5- sinqi H q(1 — cosqi) 

qi qf qf 



i + ^,COSq, 



sm qi 

qi 



(59) 
(60) 



Re 

+ 1(1 

It 



*n 



1/^0 



1 

UJ 



Im 



*n 



It 



i/^O 



1 

UJ 

8_ 

qi 



24\ 

1 2 sm qi 

qf 



COSQi 



24^ 

— ) cos qi 



(61) 



1- 



9» 



sm qt 



(62) 



Substituting Eqs. (HSJ-dM]) into Eq. ^^, one obtains for 
the real and imaginary parts of a the following expres- 
sions: 



provided that v <^ Vg- In these equations the electron 
concentration depends on the temperature as^ 



ne{T) = no 



8 {^{T)J 



(65) 



If one neglect the oscillation terms in Eq. (|64|) then we 
obtain 



sph 



i^,T) 



47rw 



-ax 



(66) 



with ujliiT) = AT:e^ne{T)/m. 



At T — > 0, Eqs. (|63)) and (|64)) agree with expressions 
known from an earlier calculations 42iii 

Another simplest calculations of conductivity tensor 
in the case of MNs with spherical shape can be made in 
approximation oi Vg ^ v with the help of the following 
formulas 



o'{uj,T)-- 
2 4 



M(r) 



4 



de 



de 



_3/2 



qe 



qt 



■ sm q^ 



(1 — COSQj 



(67) 



a"{uj,T) 



1 



1 



4e^ M(r) 

1 



de\- 



de 



.3/2 



- + — [cosq. sm qe 



(68) 



The parameter q^, (defined by Eq. ([5^ with ^ — e), is 
governed by the frequency. Depending on the ratio be- 
tween the incident frequency and the frequency i's, the 
value of g^ can be greater or less than 1. To perform sub- 
sequent analytical calculations of the integrals involved 
into Eqs. ((67)) and (l68l) . one may use the same method 
as one described above after Eq. I\'38\i . Applying it, one 
comes to the results given above by Eqs. (|55)) and[ 



B. Linewidth and "figure of merit" 



sph 



{uj 


r) = 


3e 

4" 


muj 


'2 

.q^ 


"4 

.qt 


24 

+ — 

qf 


+ 


{ 24^ 

\ qf/ 


si 



CtT 1 T 



-^ sm qi 



„ 3 e^n.jT) 

a,^^iuj,T) = -—-— 



ax 



smgj 



. sin «j ^ ( 1 — cos qi ) The linewidth of the particle-plasmon resonance is con- 

qi trolled by lifetime broadening due to various decay pro- 

1 cesses. Knowing the expressions for conductivity tensor 
J 'Eq. (|3I]) and Eq. ^^, the halfwidth of the plasmon line 
/^o-absorption can be easy calculated by means of relation 



8 
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3^ 
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1 - 




cos qi 


qi 


V 


qi) 





- + ^ cos (7i 



q^ 



smgj 



V,{uj,T) = 



4:77 



+ n^il/L,~l) 



Rea,,(w,r), (69) 



24 

72" 



cosqi 



(64) 



where n{— ^/e^) is the refractive index of the surround- 
ing media and £oo = 1 + einter is the high frequency di- 
electric constant due to interband and core transitions of 
the inner electrons in a MN's material. 
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The " figure of merit" can be expressed tlirougfi tfie (Jjj 
as well. For MNs with spheroidal shape we foundii 



nWpi(l/L/|n 
FOMnn = ^ ^^' 



1) 



(1/L 



(1) 



with geometrical factors^ 



(l + i?ll/i?i)i 



1)^(1) 



Kl 



a'-\{uj,T), (70) 



(71) 



Finally, it should be reminded that the real and imagi- 
nary parts of a dielectric constant tensor ejj can be easily 
obtained with help of tensor ajj , aa^ 



47r 



e'jJuj,T) = eoo (t'-Juj,T), 



(72) 



The high temperature of an electron gas can be 
achieved inside of MN with using laser pumping at SPR 
frequency when the hot electrons are produced. As the 
temperature of an electron gas is decreased after an ini- 
tial excitation, the SPR linewidth can increased (see 
curve 1, Fig. 1) or decreased (curve 2), depending on 
the particle radius. Such different behavior is due to the 
linewidth oscillations with the particle radius changing.— 
Earlieri2Zi2£, the results for Ag and Au nanoparticles were 
reported when the linewidth is only reduced linearly with 
the temperature lowering over a range of a crystal tem- 
peratures. 

Let us calculate the polarizability tensor for MNs at 
finite temperatures. 



IV. POLARIZABILITY TENSOR 



47r 



e''Ju;,T) = —aUu;,T). 



(73) 



The sign in the last equation and before second term 
in Eq. (j72p is defined by sign in time dependence of an 
external electric field E(i) — Eqc^*"*. The imaginary 
part of the dielectric constant will lead principally to a 
damping and broadening of SPR. 

The linewidth of SPR and FOM are strongly depen- 
dent on the particle radius. We will illustrate this depen- 
dence below building the temperature behavior of SPR 
linewidth for different R as an example. 

The SPR line width is plotted in Fig. 1 against temper- 
ature for spherical Au nanoparticles with two different 
radii, embedded in the water (n = 1.33). Calculations 
were done according to the formulas Eq. (IM)) and (p5)) 
for the following parameters of Au nanoparticlei^ vp = 
1.394xl0«cm/s, Tf = 6.41x10^ K, Wpi = 1.37x101^-1, 
and Eoo = 9.84. The electron concentration was esti- 
mated from Eq. (|65p with tiq given above by Eq. p6p . 
For Au, no ~ 5.9 x 10^^ cm^^. The behavior of T{T) is 
governed by the temperature dependence of the f" (T) in 
accordance with Eqs. (|69p and ([75)) . 




0.30 



FIG. 1. (Color online) SPR linewidth vs temperature for spherical 
Au particles with different _R (A): 50 (1) and 54.7 (2), embedded 
in water. 



The dipole electric moment for a spherical particle em- 
bedded in the media with dielectric constant e„i7 can be 
written as 



dj{Lo) = e„i 2_j '^Jfc(^) Ek{^)^ 



fc=i 



(74) 



where Ek is the k component of an external electric field 
and Ujk is the polarizability tensor of the MN. The calcu- 
lation of d is especially simple for the case, when the par- 
ticle sizes are small comparing to the some " wave length" 
^ ~ c/(-\/|e| w), which corresponds to the frequency a; in 
the particle bulk. In this case one can calculate the po- 
larizability of MN using the formulae obtained for an ex- 
ternal uniform statical field. In general terms this tensor 
can be treated as being complex value 



^i'c(w) = ajfc(w) + m^'fe(w). 



(75) 



Tensor of an electric conductivity was introduced by 
Eq. ([15]). The component of an inner electric field in 
Eq. ([T5)) is connected to a component of an external elec- 
tric in Eq. ([71]) field by expression ([3|). Because, the den- 
sity of electric current is related with the dipole moment 
by a simple relation 






(76) 



it is easy to establish the connection between polariz- 
ability and conductivity of MN. Performing the Fourier 
transformation of Eq. (|76p , we obtain 



i{^) = y d{uj)- 



(77) 



Then, using Eqs. ([74|), ([77|), and ([T5|), the polarizability 
tensor can be expressed through the conductivity tensor 



ajj by means of 



,(w,T) = - 



V 



{^,T) 



w l + Lj[ejj{uj)/e„i - 1] ' 



(78) 



Parameter V in Eq. ((78)) represents for MN with a 
spheroidal shape, the volume of spheroid {— |7ri?^i?||). 



Using Eqs. ([22]) and ^ for an isolated MN in host, 
we get the Clausius-Mossotti dipole polarizabilitj^ 



a.jj{uj,T) = 



V 



ejj{uj,T)-er^ 



AnL,e,,iLo,T) + il/L,~l)e,^ 



(79) 



where the tensor ejj would correspond to a given fre- 
quency ui, and in the case of frequencies close to the 
plasma oscillations of electrons in metal, can be presented 
as 



^Jji^^T) ~eo 



/Wpi\2 47r 



+ —iajj{uj,T). (80) 



w 



Optical absorption generally measures Ime(aj,T). 



V. LIGHT ABSORPTION AND SCATTERING 

In classical electrodynamic it is supposed that the light 
scattering results from the polarization of the scattered 
particle when it is illuminated with a beam of light. Be- 
cause the spatially uniform field causes only a dipole po- 
larization, the nanoparticle scatters the light like as a 
vibration dipole.— In the case of particles of a spheroidal 
shape, the dipole moment of the MN can be presented in 
the form^ 



d = aj_E + (a|| — a^)(m • E)m, 



(81) 



where m is the unit vector directed along the spheroid 
axis of revolution. Accounting for that the m • E = 
Ecos(p, where (p is the angle between the spheroid axis 
of revolution and the direction of an electric field E, 



d - . /^,2 c;r,2 , 



■ an cos^ (j) E. 



(82) 



Using our forerunning formulae for d and a, the tem- 
perature dependence of light absorption and scattering 
can be calculated. The overall absorption and scattering 
cross-sections for a spheroidal MN can presented as^ 



Cabs = 47r — ^e^ [(ImQ;||)cos^0+ (ImQ;j^)sin^ 



(83) 



Csca^-^Tri^-V^i) (laiipcos^^+laj^psin^^). (84) 

The square of diagonal components of the polarizability 
tensor can be presented in the Lorenzian form as 



V 



\iy 



2 m - e,„)cj2 - ^2^ J + (2W7(|-) 



and the imaginary part of the polarizability tensor is 

,3i 



Ima 



V 



(1) - 4^ 



2a;-^^,„7/|N 



(")/ ftj2_cj2 1 +i^2ujj^ii 



Here 



Em + Lni\{^ — Em) ' 



11) 



^?M = Wni^(^) 



11) 



(86) 
(87) 

^liiT) (88) 



is the II and _L surface plasmon frequencies, and 

2TTLni\ 
7(1) ^ 7(1) (-,T) = --_^_i__^Rea(,)(a;,T) 

(89) 
is the half- width of the resonance for light polarized along 
(II) or across (_L) the rotation axis of the spheroid, cr||, 
a± are the corresponding components of the conductiv- 
ity tensor, and Ly, L± are the geometrical factors given 
above by Eqs. (|7T]) . 

On resonance, the polarizability tensor looks as 



a(|)(a;,r) 



V 



£,m^ni\{T) 



8ttL 



(|)7(|)(^,T) 



(90) 



7/ II \ iaK.eii ai uj ^ ujn\\. 

resonant frequencies is so small that their input into the 
scattering cross section can be neglected. 

In the case of MN with a spherical shape Ly — L± — 
1/3, and Eqs. (IH71)-(IHni) become the forms 



e« 



3e„ 



2e™ + 1 



W(|)(r) = w. 



^Ut) 



sp 



7(|)(a;,r) = 7,.p(a;,r) 



2e„ 



2tt 



(91) 



(92) 



2e„,, + 1 



a{io,T). (93) 



Then the expressions ((84)) and ((83)) can be easily trans- 
formed to the forms 



Cabs{L^,T)^Wefl^- 



e"iu:,T) 



AttL 



(1). 



(1)) +(2-7(1) 



c[e'iuj,T) + 2emf + [e"iu;,Tf' 

(94) 

V ( T^-3y2 ^ .^^4 [,>,T)-e„]2 + [,"(a.,r)]2 

^sca[^7 J ) — -T, ^m 1 " 



(85) 



27T ^'^\cJ [e'{Lo,T) + 2e„,f + [e"{u,T)f'' 

(95) 
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where we have combined Eqs. (|9H) - ((93)) with the real 
and imaginary parts of Eq. ^fSU\i . The most intense cross- 
sections are observed at the frequency which corresponds 
to the plasmon resonance of a spherical nanoparticle in 
the vacuum. Formulae (|95|) and Eq. (p4|) yield a reso- 
nance when e'{uj,T) = — 2em- The absorption and the 
scattering cross-section defined by Eqs. ([55)) and ^E^ 
reach the maximal values, as well, at the frequencies of 
the SPRs uj — LOf \\\. Near the surface plasmon resonance, 

light may interact with the particle over a cross-sectional 
area usually (at low temperatures) larger than the geo- 
metric cross section of the particle because the polariz- 
ability of the particle becomes very high in this frequency 
range. 

We will calculate the the scattering efficiency of the 
MN, which we define as the ratio 



C 



'e/ 



a 



Ca 



geom 



C 



(96) 



where C„r„„,, is the geometrical cross-section of an indi- 



^ geom 



vidual particle and C is defined by Eqs. 
by Eqs. ([M ]) and (|95l) . 



and 
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FIG. 2. (Color online) The absorption crossection vs temperature 
for spherical Au particles with different R (A): 78 (1), 100 (2), 125 
(3), and 150 (4), embedded in water. 

Figure 2 shows the variation of an absorption efficiency 
as a function of temperature for a gold nanoparticles with 
a spherical shape embedded in the water (n — 1.33). 
The calculations were performed in an accordance with 
Eqs. (|96)) and (|94|. with the use of the same parame- 
ters as above for Fig. 1. For zero temperature we get 
the classical result,-'- when the SPR intensity depends 
on the particle radius: larger particles have a larger scat- 
tering or absorption cross sections. The situation is dras- 
tically changed at finite electron temperatures. As can 
see from Figure 2, the absorption efficiency can not only 
increased (right hand sides of curves 1-4) with tempera- 
ture decreasing, but decreases as well (left hand sides of 
curves 1-4). We have an interesting situation, when the 
absorption crossection with the temperature lowering is 
increased at first, reaches its peak, then is decreased, and 



start to increase again. These observations are explained 
with the conductivity of MN, which oscillate with parti- 
cle radius changing.— 

In the case of light scattering (see Fig. 3), the effect 
of MN size is more pronounced because the scattering 
cross-section is proportional to the square of a particle 
volume. 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 

FIG. 3. (Color online) The same as in Fig. 2, but for scattering 
crossection. 

From Fig. 3, one can see that the efficiency of light 
scattering at resonance frequency has not only increased 
(curves 1-3) with the enhancement of the MN radius, 
but can be decreased (curve 4) as well. The latter means 
that the efficiency of scattering at resonance frequency 
can be suppressed with temperature rise for some par- 
ticle radiuses. Depending on the electron temperature, 
the scattering efficiency reaches the maximum value, the 
peak position of which depends on the nanoparticle ra- 
dius. 

The peak positions of the resonance plasmon absorp- 
tion (or scattering) by Au nanoparticle embedded in wa- 
ter can shifted toward the shorter wavelengths as the 
temperature of electron gas is lowered after an initial 
excitation. This behavior is consistent with the known 
temperature-dependent shift of the differential transmis- 
sion peak in isolated sodium nanoparticle j^ 



VI. DISCUSSION OF RESULTS 

Here we discuss shortly the effects caused by devia- 
tions from the spherical particle shape. When the sym- 
metry of the particle decreases, the number of resonance 
peaks increases. If the shape of particle deviates from 
spherical one, the 3-fold degeneracy dipole mode split 
into two (for spheroid) or three (for ellipsoid) modes giv- 
ing rise to correspondent scattering peaks. In this case 
the plasmon response depends strongly on the particle 
position relative to the illumination direction. To eval- 
uate the electron temperature effect on the absorption 
or scattering crossections for various illumination direc- 
tions, one can perform a numerical calculations with the 
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use of Eqs. (|83l)-(189l) and Eq. ((3T|) . 

The particle absorption or scattering effects depend on 
the size and shape of the individual particles. Metallic 
particles that are much smaller than the wavelength of 
light tend to absorb more and hence absorption is dom- 
inated. The scattered light intensity from small MNs is 
extremely low. Indeed, Rayleigh scattering theory pre- 
dicts that the scattering cross section of a spherical clus- 
ter should drop as the sixth power of the diameter. 

The quantity that is measured experimentally is 
frequency-resolved time-dependent transmission^"^ To„ 
after the exitation of the MN by a pump laser pulse. 
Then, the static equilibrium transmission To// of the 
sample measured in the absence of the pump is sub- 
tracted to obtain the relative differential transmission^^ 
(5T/T ={Ton — '^off)/^off- In the thermal equilibrium, 
the 5T can be related to the temperature dependent ab- 
sorption cross section by^i 



^T 
T 



27ri?2 



[Cabs{^,T) - Cabs{<^,Ta)\, 



(97) 



where Tq is the temperature of of the MN before excita- 
tion. This relation holds for very small sizes of MN when 
the reflectivity of the MN is very small. 

A direct comparison of theoretical results with most of 
the available experimental measurements of the optical 
properties of MNs is still a matter of debate because in- 
homogeneities in nanoparticle size, shape, and local en- 
vironment hide the homogeneous width of the surface 
plasmon resonance. 



VII. CONCLUSIONS 



close to surface plasmon frequencies of a spheroidal MN 
is studied in detail. 

The high sensitivity of the temperature dependence of 
SPR linewidth to the radius of a particle was established. 
Even small variation in the particle radius can drastically 
change the run of the SPR linewidth curve vs tempera- 
ture. This is due to the linewidth oscillations when the 
radius of MN is varied. The efficiencies of light absorp- 
tion and scattering for various temperatures strongly de- 
pend on the particle radius as well. The An nanoparticles 
with different radii are used for illustration. 
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Appendix A 

Below we advance the values for the integrals appeared 
in calculations of Eq. (P^ . Sec. Ill: 



-1 
1 /7 

2\^y V 2Va J "" V 2^/a 



Erfrii±I^VErf^^-2^ 



(Al) 



We have applied a kinetic theory to calculate the op- 
tical properties of a metal nanoparticle embedded in a 
dielectric media at different temperatures. It allows to 
determine the cross-section of light absorption or scatter- 
ing for various polarizations of the incident electromag- 
netic wave. The obtained analytical formulas provides to 
evaluate the dynamic of a light absorption or scattering 
intensities at the plasmon frequencies for nanoparticles 
with different radii and shapes, when the temperature 
of electron gas becomes settle after an action of a laser 
pulse. 

The analytical expressions for electroconductivity and 
polarizability tensors have been obtained for the case 
of MNs with a spherical shape. They allow to calcu- 
late the number other physical quantities (e.g., reso- 
nance linewidth, figure of merit etc) for a given temper- 
ature. We studied theoretically the variation of both the 
linewidth of SPR and the efficiency of the light absorption 
and scattering by MN with changing of the electron tem- 
perature. The case where the frequency of a laser beam is 



dz z^e-('+^y = ^(1 + 2a) e-^% (A2) 



provided that —tt/A < a < n/A. 



oo 

/ dz z^ Eric(z + - 



a\ 1 / 3 , _4„ 



provided that Re a > 0. 



(A3) 



OO 



= - I a + a 



(A4) 
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